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A METHOD, AN INSTALLATION, AND A COMPUTER PROGRAM FOR 
ESTIMATING THE INITIAL SIZE OF A POPULATION OF NUCLEIC 
ACIDS, IN PARTICULAR BY PCR 

The present invention relates to estimating the 
initial size of a population of interest, in a sample 
subjected to a succession of amplification reactions. 

BACKGROUND OF THE INVENTION 

The present invention finds a particularly 
advantageous, but non-limiting, application in 
determining an initial quantity of nucleic acids in a 
sample subjected to a polymerase chain reaction (PCR) in 
real time. A technique of this type, known as 11 PCR 
quantification 11 , is used in particular for evaluating the 
number of copies of pathogenic agents (e.g. of the human 
immunodeficiency virus (HIV) ) in a sample of body fluids 
taken from a patient, typically in the context of a 
medical checkup. 

Reference is made to Figure 1 for a brief 
description of the diagrammatic appearance of a real time 
PCR amplification curve with PCR cycle index numbers 
plotted along the abscissa and, in the example shown, 
with quantities of fluorescence emitted (in arbitrary 
units) as measured for each PCR cycle plotted up the 
ordinate. For each PCR cycle, it should be understood 
that the sample is subjected to temperature variations 
enabling DNA polymerase to amplify nucleic acids and 
enabling the corresponding PCR products to be detected by 
fluorescent molecules. By plotting the measured 

fluorescence F n as a function of PCR cycle number n, 
variation is obtained of the type shown in Figure 1, and 
comprises at least: 

• a first portion BN where fluorescence measurements 
coincide substantially with the background noise of the 
apparatus for measuring fluorescence; 

a second portion EXP in which the measured 
quantities of fluorescence increase in substantially 
exponential manner; 



• a third portion LIN in which the increase in the 
measured quantities of fluorescence is significantly 
attenuated and behaves overall in substantially linear 
manner; and 

a fourth portion PLA in which fluorescence 
measurements reach a plateau stage. 

It should be observed that for the initial PCR 
cycles (first and second portions), the population of 
interest increases in substantially exponential manner, 
whereas for the following cycles (third and fourth 
portions) , other phenomena come into competition with 
growth in the population of interest, so that said growth 
is then damped up to the plateau stage PLA. 

The document "Mathematics of quantitative kinetic 
PCR and the application of standard curves" by 
R.G. Rutledge and C. Cote, published in Nucleic Acids 
Research, 2003, Vol. 31, No. 16, discloses a method of 
estimating the unknown initial quantity of nucleic acids 
in a sample of interest by means of PCR. That method 
consists in using a plurality of samples having known 
initial quantities of nucleic acids, referred to as 
"standards", in order to determine by interpolation the 
initial quantity of nucleic acids present in the sample 
of interest . 

In general, the greater the initial quantity of 
nucleic acids in a sample, the sooner a detectable 
quantity of PCR product is obtained, i.e. the sooner a 
detectable quantity of emitted fluorescence is obtained. 
With reference to Figure 2, relating to the prior art, it 
will be understood that the initial population in the 
standard Stl is greater than that in the standard St2 
which is greater than that in the standard St3, etc., 
since the cycle Ctl for the standard Stl occurs before 
the corresponding cycle Ct2 for the standard St2, which 
occurs before the cycle Ct3 for the standard St3, etc. 

Thus, such a CT cycle, corresponding to the cycle at 
which the fluorescence measurements reach a fluorescence 
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threshold THR (as shown in Figure 2), sets at an 
arbitrary level (typically below the background noise) , 
and acts as a parameter representative of the initial 
size N 0 of a population of nucleic acids subjected to the 
5 PCR cycles. Use has been made of this observation in the 
above-cited prior art to establish a relationship of the 
kind shown in Figure 3 between cycle numbers Ctl, Ct2, 
Ct3, Ct4 for a plurality of standards having known 
initial populations, and their initial populations N 0 x f 

10 N 0 2 , N 0 3 , N 0 4 . Thus, by plotting the cycles Ctl, Ct2, Ct3, 
Ct4, etc. up the ordinates and the logarithm of the 
initial population sizes Nq 1 , N 0 2 , N 0 3 , N 0 4 along the 
abscissa, a regression slope REG is obtained. On this 
regression slope PEG, the cycle Ctint detected for the 

15 sample of interest- is plotted (dashed-line arrow Fl) . By 
interpolation on the regression slope REG (dashed-line 
arrow F2) , the initial population size N 0 int is then 
determined for the sample of interest. 

Although that method is in widespread use, it 

20 nevertheless presents some drawbacks. 

Firstly, it requires the use of a plurality of 
standard samples having respective known initial 
populations . 

Secondly, the method depends on the judgment of the 
25 user, since the fluorescence threshold value, as selected 
by the user, has a direct influence on the values of the 
Ct cycles in the amplification curves, and consequently 
on the estimated values for the initial population size 
in the sample of interest. The threshold value also has 
30 an impact on the accuracy of the result, since accuracy 
is generally better if the threshold is selected to lie 
in the exponential growth stage EXP of the amplification 
curve. Nevertheless, in practice, it is difficult for 
the user to know whether the fluorescence threshold level 
35 THR that has been set does indeed correspond to the 
exponential stage of the curves, and does so for all of 
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the samples (the standard samples and the sample of 
interest ) . 

Finally, the method assumes without any verification 
that the population has the same amplification yield in 
5 the sample of interest and in all of the standard 
samples. Thus, if the sample of interest contains PCR 
inhibitors, as is typically the case, then its result 
will be falsely lowered. 

It should thus be understood that the prior art 

10 technique depends on the fluorescence threshold THR as 
defined by the user. The value selected has an influence 
on the values of the Ct cycles and consequently on 
determining the initial quantity in the sample of 
interest. That is one of the reasons why a large amount 

15 of work has recently been undertaken to automate Ct cycle 
detection and make it reliable. 
OBJECTS AND SUMMARY OF THE INVENTION 

The present invention seeks to improve the situation 
by proposing an approach that is completely different. 

20 Firstly, the invention provides a method, the method 

being implemented by computer means for quantifying in 
absolute and/or relative manner an initial population of 
nucleic acids in a sample of interest. The sample is 
subjected to a succession of applications of a reaction 

25 for amplifying the population of interest. In very 
general manner, this amplification may be undertaken by 
implementing successive PCR cycles, however any other 
amplification technique could also be used. Above all, 
it should be understood that the amplification needs 

30 merely to be defined by a reaction yield, as described 
below. During these successive amplification operations, 
experimental measurements are taken that are 
representative of a current population size, at least in 
the sample of interest. It will be understood that one 

35 or more measurements can be taken after or during each 
amplification reaction without loss of generality. 
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In a presently preferred definition of the 
invention, the method in the meaning of the invention 
comprises the following steps: 

a) providing a model of the yield of the 
5 amplification reaction as a function of the succession of 

amplifications, said model comprising: 

a substantially constant stage for a first 
portion of the applications of the amplification 
reaction; and 

10 • a non-constant stage for a second portion of 

the applications of the amplification reaction; 

the first and second portions being united by a 
changeover region in which yield changes over between the 
constant and non-constant stages, said region having an 

15 amplification index corresponding substantially to the 
changeover ; 

b) using the yield model to express a relationship 
involving at least the changeover index and a parameter 
representative of the initial population size in the 

20 sample of interest; 

c) determining at least the changeover index by 
comparison with the experimental measurements; and, in a 
subsequent or immediately following step d) deducing 
therefrom the initial population size in the sample of 

25 interest. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Other advantages and characteristics of the 
invention appear on reading the following detailed 
description of an implementation given below by way of 
30 example with reference to the accompanying figures, in 
which : 

• Figure 1 relates to the prior art and represents 
variation in the measured quantity of fluorescence as a 
function of the number of PCR cycles, as described above; 
35 • Figure 2 relates to the prior art and is 

representative of the increasing quantities of 
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fluorescence that are emitted as a function of the number 
of PCR cycles, as described above; 

■ Figure 3 represents an interpolation method for 
determining the initial quantity of the population of 
5 interest in the sample of interest using a method known 
in the prior art, and described above; 

Figure 4A is a diagram showing variation in the 
above-described experimental measurements as a function 
of the succession of amplifications applied to the sample 
10 of interests- 
Figure 4B is a diagram showing variation in the 
yield of the amplification reaction, obtained from 
experimental measurements, as a function of the 
succession of amplifications applied to the sample of 
15 interests- 
Figure 5 plots a regression relationship between 
the yield changeover indices that occur at the changeover 
between the constant stage and the non-constant stage, 
and the logarithms of the initial populations for 
20 standard samples and for the sample of interest, for use 
in a first implementation; 

• Figure 6A shows typical variation in the measured 
quantity of fluorescence after it has been adjusted by 
taking account of background noise specific to the 

25 measurements, and plotted as a function of the number n 
of PCR reaction cycles; 

• Figure 6B shows the variation in the effectiveness 
of the PCR shown in Figure 6A as a function of the number 
n of cycles; 

30 • Figure 7 compares the experimental variation in 

emitted fluorescence as shown in Figure 6A with the 
results obtained by applying the emitted fluorescence 
model obtained by including an effectiveness model in a 
second implementation; 

35 ■ Figure 8 is a comparison between the variation in 

the effectiveness of Figure 6B and the application of an 



effectiveness model deduced from the emitted fluorescence 
model of Figure 7; 

• Figure 9 is a flow chart outlining the main steps 
in the method in a proposed implementation of the present 
invention; and 

Figure 10 is a diagram of an installation for 
quantifying the initial population of a sample of 
interest . 

MORE DETAILED DESCRIPTION 

Reference is made to Figures 4A and 4B for briefly 
describing a few principles of the invention illustrating 
the characteristics of the above method. 

Firstly, it should be understood that Figure 4A 
plots a succession of experimental measurements F n 
representative of - the current size of a population of 
interest which is being subjected progressively to a 
succession of amplification reactions, each reaction 
being indexed by an index number n. In the non-limiting 
example described herein, this succession of reactions 
corresponds to a succession of PCR cycles. In this non- 
limiting example, the experimental measurements F n 
correspond to measured quantities of fluorescence on each 
PCR cycle. Thus, in a quantification method which 

combines the PCR reaction and the fluorescence emitted by 
the sample of interest, fluorescent reagents are 
introduced into the sample so that the fluorescence that 
is emitted during a PCR cycle is proportional to the size 
of the nucleic acid population in the sample. Indeed, it 
can be preferable to perform a plurality of measurements 
or no measurements at all for certain PCR cycles. 
Furthermore, more generally, the measurement method may 
make use of techniques other than fluorescence, even if 
fluorescence is the method that is often used for 
quantifying by PCR. Finally, it should be understood 
that other amplification techniques could be implemented 
in the context of the present invention, providing it is 
possible to track variation in the yield of the reaction 
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corresponding to the amplification. Since the example 
described below relates preferentially to PCR cycles, 
reference is made to the "effectiveness of the PCR" 
written E n for each PCR cycle of index n, in order to 
5 refer to the yield of the amplification reaction. 

As mentioned above with reference to Figure 1, 
Figure 4A mainly comprises two regions in which: 

• during the initial PCR cycles (portion EXP) , the 
population increases substantially exponentially; whereas 

10 * during the following cycles (the LIN and PLA 

portions), other phenomena come into competition with 

growth of the population of interest, so the growth 

becomes damped. 

The following two assumptions are made: 
15 • the yield of the reaction E n is relatively constant 

during the initial cycles over the portion EXP; and 

after some number of cycles have been performed, 

the yield E n of the reaction starts to decrease over the 

portions LIN and PLA. 
20 This decrease in yield may have a variety of 

explanations, in particular a degradation and/or a lack 

of PCR reagents (DNA polymerase, dNTPs, primers, etc.) 

and/or inhibition by the products that are made 

themselves . 

25 It is assumed herein that the yield is initially 

constant and that it subsequently decreases. 
Nevertheless, it should be understood that the invention 
applies more generally to the context of yield: 

• that is initially constant, which corresponds to a 
30 normal situation for growth by amplification; and 

that is subsequently not constant (decreases or 
increases) which corresponds to a situation that is 
substantially abnormal . 

In the context of reactions for amplifying the 
35 quantity of nucleic acids, it has been found that the 
yield often changes over from a constant stage to a non- 
constant stage. In the meaning of the invention, 
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advantage is taken of this observation to deduce 
therefrom the initial quantity of nucleic acids, as 
described below in detail. Initially, it is merely 
stated that the yield can also change over from a non- 
5 constant stage during early cycles to a subsequent 
constant stage. The present invention is equally 

applicable to such a circumstance. In general, it should 
therefore be understood that in the meaning of the 
invention, a changeover of yield between a constant stage 

10 and a non-constant stage is detected. 

The objective is to find the initial size of the 
population that has been subjected to amplification. 
With reference to Figure 4, it will be understood that 
the measurement F 0 representative of this initial 

15 population size, -which coincides in practice with the 
measurement background noise BN, cannot be used on its 
own for determining directly the initial population size. 
In the prior art, attempts have been made to quantify 
this initial population size by making use of the 

20 exponential stage, i.e. a stage that occurs typically on 
exiting background noise. A threshold cycle Ct is then 
determined (corresponding to point PA for "prior art") in 
Figure 4A. As mentioned above, in this region 

measurements are often affected by noise and it is 

25 difficult to determine accurately a threshold cycle Ct 
representative of exiting background noise. 

In a completely different approach, the present 
invention instead makes use of nearly all of the points 
of the amplification curve in order to determine 

30 accurately a region CHO where the yield changes over 
between a constant stage and a non-constant stage, 
typically in present circumstances between the 
exponential stage EXP and the linear stage LIN. It will 
be understood that measurements are logically less 

35 affected by noise in this region CHO than in the 
background noise exit region since the region CHO occurs 
during later cycles. Furthermore, particularly because 
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of the mathematical properties associated with yield, it 
is shown below that, most advantageously, the number of 
standards that need to be used for quantifying the 
initial size of the population of interest is smaller 
than the number of standards used in prior art 
quantification . 

The relationship for associating the changeover 
region CHO with the initial size of the population of 
interest is briefly described below. The yield of an 
amplification reaction is given by: 
N n+ i = N n + E n x N n 
in which: 

• N n is the size of the population of interest after 
an amplification of index n in a succession of 
amplifications; 

N n+1 is the size of the population of interest 
after a following amplification, of index n+1, in the 
above-mentioned succession of amplifications; and 

• E n is the yield of the amplification reaction of 
index n in the above-mentioned succession of 
amplifications . 

Reformulating this relationship as a recurrence 
relationship, we obtain: 

N n+1 = (1 + E n ) (1 + E n -i) (1 + E„- 2 ) ... (1 + E 0 )N 0 

where N 0 is the initial size of the population of 

interest. So long as the yield E n is constant, it will be 

understood that the above relationship can be written 

more simply as follows: 

N n+ i = N 0 x (1 + E 0 )n+i 

where the index n+1 has not yet reached the changeover 
region CHO. While the yield is constant during the 
initial cycles, the following applies: 

En = En-l = E n _ 2 = ... = Eo 

where E 0 is the value of the yield during the constant 
stage. Nevertheless, when the index n+1 moves into the 
changeover region CHO, the relationship becomes: 
N n +i = N 0 x (1 + E 0 ) C eep x function (C EE p, n+1) 
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where : 

(Ceep ~ 1) ^ s the last index of the amplification 
reaction during which the yield is still constant (it 
will thus be understood that the index C EEP itself 
5 represents the changeover index proper between the 
exponential stage and the linear stage) / and 

the term function (C EEP , n+1) is a particular 
function characterizing the non-constant stage of the 
yield and that depends at least on the changeover index 

10 C EEP and on the current amplification index n+1. 

It can thus be seen how it is possible to associate 
the changeover index C EEP and the initial size N 0 of the 
population of interest. At this stage it can be 

understood that steps a) and b) of the above-defined 

15 method have already been implemented. 

A first implementation consists in determining the 
changeover index C EEP experimentally and in correlating it. 
with the initial size by regression by using a plurality 
of standard samples that are subjected to the same 

20 amplification treatment as the sample of interest. Under 
such circumstances, it will be understood that steps b) 
and c) of the above-defined method are merely 
interchanged since initially the changeover index C EEP 
(step c) ) is determined experimentally, and subsequently 

25 the relationship between the index C EEP and the initial 
size N 0 (step b) ) is determined in order to end up with 
the initial size N 0 (step d) ) . 

Before describing all of these steps in detail in 
the meaning of the first implementation, a method is 

30 described for determining the index C EEP on the basis of 
experimental measurements. In particular, it will be 
understood that this method of determining the index C EEP 
experimentally can be applied to another implementation 
that is different from the above-mentioned first 

35 implementation . 

Returning to the relationship between the 
effectiveness E n of a given cycle n and the current size 
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of the population of interest in the same cycle N n and in 
a subsequent cycle N n+1 , the effectiveness of the 
amplification can be expressed as follows: 

En = < N n + l/N n ) - 1 

5 In certain circumstances, in particular when there 

is no need to take account of background noise BN in the 
measurements, it is possible to a first approximation to 
assume that the measurements are substantially 
proportional to the current size of the population of 

10 interest. Nevertheless, in practice, account will more 
often be taken of measurement drift, with corrected 
experimental measurements F' n being determined on the 
basis of direct measurements F n as shown in Figure 4A. 

A prior step of processing the experimental 

15 measurements F n -is preferably applied, this step 
consisting in subtracting the background noise BN and 
subsequently in introducing compensation to take account 
of a non-zero measurement 8 representative of the initial 
population size. In the example shown in Figure 4A, the 

20 variation in the background noise BN as a function of the 
index n can be represented by a linear function since 
tests have shown that a linear model is satisfactory for 
fluorescence measurements in PCR. Nevertheless, in 

certain circumstances it may be preferable to use an 

25 exponentially-varying model. In any event, a model is 
applied that complies best with variation in the 
background noise BN as given typically by the initial 
measurement points. Thereafter, the selected model for 
variation in background noise BN is subtracted from all 

30 of the experimental measurement values F n . By applying 
this step, it will be understood that the theoretical 
fluorescence measurement F 0 is reduced to a measurement 
value of zero, corresponding to an initial population 
size N 0 of zero, which is not representative of physical 

35 reality. Consequently, it is advantageous to apply 

compensation for this correction as follows: 
F 1 n = F n - BN + 8 
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where : 

• the term F f n corresponds to a corrected measurement 
for a current index n; 

the term F n corresponds to the raw experimental 
5 measurement at said current index n; 

the term BN corresponds to the value for the 
background noise as modeled for the current index n; and 

• 8 is the corresponding compensation term which is 
assumed to be constant in the example being described and 

10 which directly represents the initial population size N 0 . 

Although these steps of correcting for background 
noise are very advantageous in determining the changeover 
index C EEP , they may also be applied to any determination 
and quantification of the initial population size N 0 

15 whenever background noise is likely to falsify 
measurement of said population size N 0 . In this respect, 
these steps may constitute the subject matter of separate 
protection, where appropriate. 

The corrected measurements F f n as obtained in this 

20 way are advantageously proportional to the current 
population sizes N n in the samples of interest, such that 
the yield E n can now be expressed directly as a function 
of measurement values (corrected as described above) , by 
the following relationship: 

25 E n = (F" n+1 /F' n ) - 1 

Thus, from the experimental measurements F n of 
Figure 4A, corrected measurements F' n are obtained from 
which there is subsequently determined the variation in 
the effectiveness E n as a function of the succession of 

30 indices n, as shown in Figure 4B. 

In short, the experimental measurements are 
expressed in the form of an experimental variation in the 
effectiveness E n of the kind shown in Figure 4B as a 
function of the succession of amplifications n. This 

35 provides an experimentally-determined variation for the 
yield comprising: 
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a perceptibly noisy first region for low 
amplification indices n (specifically prior to the cycle 
CG in the example of Figure 4B) ; and 

• followed by a second region exhibiting less noise 
5 for higher amplification indices (at least after the 
changeover region CHO) . 

At least in the most usual circumstance of 
amplification by PCR and measurement by fluorescence, the 
non-constant stage of yield is decreasing and corresponds 

10 to said second region presenting little noise (as shown 
in Figure 4B) . Specifically for the purpose of 

eliminating measurement points that run the risk of 
falsifying results when selecting a model to apply to the 
variation in yield: 

15 • a crude value E 0 is estimated for the constant 

yield stage; and 

particularly when searching for the changeover 
index C EEP , at least some of the measurements in the less 
noisy second region are ignored for which the estimated 

20 yield is less than a threshold value, e.g. corresponding 
to some fraction of the constant stage E 0 . 

These points NEG (Figure 4) that are eliminated are 
typically those that correspond to very high 
amplification indices n and that might no longer satisfy 

25 the model for effectiveness which is selected 
substantially around the changeover region CHO. By way 
of example, in order to eliminate them, an average is 
evaluated for the constant stage of yield E 0 , typically 
for the initial indices n. Thereafter, a threshold value 

30 is selected that corresponds to a fraction of the average 
found for the constant stage E 0 , e.g. 10%. Thereafter, 
starting from the highest indices n, all measurement 
points NEG of measured yield lower than or equal to said 
threshold value are eliminated. This step, which is most 

35 advantageous for detecting the index C EEP , can 
nevertheless be applied to any determination based on 
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yield E n , and may constitute the subject matter of 
separate protection, where appropriate. 

When yield presents a non-constant stage in which 
yield is decreasing and which follows a constant stage, 
5 as shown in Figure 4B, the changeover region CHO is 
identified by working in the direction of decreasing 
index numbers n, starting from the less noisy second 
region, and by detecting a coarse index CG for which the 
yield passes a predetermined value. Thus, with reference 

10 to Figure 4B, going in the direction of decreasing index 
numbers n so as to rise towards the changeover region 
CHO, the yield associated with each index number is 
evaluated. For the first measurement point of yield that 
is significantly greater than the above-mentioned 

15 predetermined value, it is considered that the above- 
mentioned coarse index CG has been detected and 
corresponds to the index of the measurement point. 

As described below in a subsequent implementation, 
it is possible for each measurement point to model the 

20 variation in its yield as though said set point itself 
corresponds to the changeover index C EEP . In that 

implementation, if the constant yield stage E 0 is 
estimated, and if the estimated value exceeds the above- 
mentioned predetermined value, then the point is 

25 considered as corresponding to the coarse index CG. 

In general, a maximum yield has a value of 1 so it 
is possible to select the above-mentioned predetermined 
value as being equal to 1. Nevertheless, this can be 
varied, and, for example, provision can be made to set 

30 the predetermined value as corresponding to the mean 
yield E 0 as evaluated over the initial reaction cycles. 

Thereafter, the estimate of the value of the 
amplification index C EEP in the changeover region is 
refined, which value may advantageously be a fraction, by 

35 working in the direction of increasing amplification 
index numbers, starting from the coarse index CG, and by 
detecting an amplification index for which the yield is 
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approximately equal to the above-mentioned predetermined 
value. Thus, referring again to Figure 4B, in order to 
refine the search for the changeover index C EEP after the 
coarse index CG has been determined, a search is made 
5 downwards starting from the coarse index CG and going in 
the direction of increasing index number n, in steps of a 
size smaller than one whole index, and the abscissa value 
is determined, e.g. by interpolation, at which the 
predetermined value is crossed. Typically, so long as 

10 the constant stage value remains greater than 1, the 
search is continued in the direction of increasing number 
n, and the index (C EEP - 1) preceding the changeover is 
determined as soon as the constant value E 0 is equal to or 
very close to 1. That is why it is appropriate to select 

15 a search step si?e corresponding to a fraction of the 
index, for example 10% of one cycle n. 

In the above-mentioned first implementation, a 
plurality of standard samples are provided having 
respective known initial population sizes, and the 

20 succession of amplifications is applied thereto under 
substantially the same conditions as for the sample of 
interest. Their respective changeover indices are 

determined in accordance with above-described steps a), 
b) , and c) . In step d) : 

25 • a dependency relationship is established between 

the initial population sizes of the standard samples N 0 st 
and their indices C EEP st ; and 

• after determining the index C EEP for the sample of 
interest, the initial size of the population of interest 

30 N 0 is determined by interpolation on that dependency 
relationship . 

Thus, with reference to Figure 5, a dependency 
relationship can be established between the changeover 
cycles C EEP l , C EEP 2 of the standards and their initial 
35 concentrations Nq 1 , N 0 2 (actually the logarithms thereof) , 
e.g. by regression. By measuring the changeover index 
C EEP int for the sample of interest and by plotting its 
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value on the regression slope of Figure 5, the initial 

concentration N 0 int in the sample of interest is obtained 

by interpolation. 

This first implementation is thus quite similar to 
5 that of the prior art described with reference to 

Figure 3. Nevertheless, it should not be forgotten that 

the changeover index C EEP on which this first 

implementation relies does not correspond in any way to 

the threshold cycle Ct of the prior art. 
10 In an approach that is significantly different from 

this first implementation: 

■ in step b) , use is made of the yield model to 

express variation that is parameterized as a function of 

the succession of amplifications, said variation making 
15 use of at least one parameter representing the changeover 

index C EEP ; and 

• in step c) , at least said parameter representing 

the changeover index C EEP is determined by comparison with 

the experimental measurements. 
20 In a second implementation, this parameterized 

variation is representative of the current population 

size N n in the sample of interest. 

Typically, this parameterized variation can be drawn 

from an expression of the type given above: 
25 N n+ i = N 0 x (1 + E 0 ) C eep x function (C EEP , n+1) 

Thus, in addition to a parameter representing the 

changeover index C EEP , this variation makes use of a 

parameter representative of the initial population size N 0 

in the sample of interest. 
30 Thereafter, in steps c) and d) of this second 

implementation, these two parameters C EEP and N 0 are 

determined substantially together. 

Previously, in step a) , it is necessary to determine 

a model for the above-mentioned function function (C EEP , 
35 n+1) . 

Usually, for PCR quantification, a model is selected 
for the non-constant stage of the yield corresponding to 
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a decreasing exponential having a decrease parameter p 
which is described in greater detail below. This 
decrease parameter p is then determined in step c) , at 
least with the changeover index C EEP , by comparison with 
the experimental measurements. 

Thus, in this second implementation, once the yield 
model E n has been selected, it is applied to the general 
expression for the current population size N n given by the 
above relationship. This provides a model for variation 
in the current population size N n . 

Nevertheless, unless the experimental measurements 
give the value for the current population size N n directly 
(which is rarely true in practice at present), it is 
appropriate subsequently to model the experimental 
measurements F n -themselves, taking account of the 
subtracted background noise and the subsequent 
compensation e as described above. 

Thus, in a presently preferred implementation, the 
above-mentioned parameterized variation: 

is representative of experimental measurements; 

and 

includes a parameter corresponding to a 
measurement value F 0 representative of the initial 
population size. 

Thereafter, the measured value of the initial 
population size F 0 is determined by comparing said 
parameterized variation F n with the experimental 
measurements . 

In order to perform this comparison, it is possible, 
for example, to adjust the parameters F 0 , E 0 , C EEP and the 
decrease parameter p in the model of the measurements F n 
by using statistical correlations (typically the least 
squares method) applied to the raw experimental 
measurements. An example implementation is described in 
detail below. 

Initially, variation is obtained for the measured 
and adjusted quantity of fluorescence as a function of 
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the number of PCR cycles that have been applied, as shown 
for example- in Figure 6A. This figure shows the 

amplification curve for a sample of interest containing 
nucleic acids, in this case a fragment of DNA having an 
5 initial quantity of 100,000 copies, marked by the 
SYBRGREEN intercalant during the PCR reaction which is 
performed on the I-Cycler IQ® apparatus from the supplier 
BI-RAD® . 

In the example described, it will be understood that 
10 the amplification reaction is a PCR reaction in real 
time. The experimental measurement represents quantities 
of emitted fluorescence. 

The fluorescence of cycle n after adjustment for 
background noise, as described above, is written F n below. 
15 The theoretical initial fluorescence before the first 
cycle is written F 0 . The effectiveness of the PCR in 
cycle n is written E n . The total number of cycles 
performed during the PCR reaction is written N . 

By assumption, the fluorescence measured on each 
20 cycle n of the PCR reaction cycle is defined by: 

F n + i * F n (l + E n ) for all n e {0 , 1, 2,...,N -l} (1) 
with 0 < E n < 1. 

The effectiveness of the reaction on each cycle n is 
calculated as follows: 

25 E n =^- 1*^-1 for all n e {0 , 1 , 2 ,...,N - 1} (2) 

It should be observed that equation (1) is assumed 
to be true for n=0 . Nevertheless, by definition, the 
initial fluorescence F 0 is unknown. It is therefore not 
possible to calculate the effectiveness on the first 

30 cycle E 0 directly for formula (2) . 

Figure 6B shows the effectiveness of the PCR 
reaction as approximated by formula (2) and on the basis 
of the adjusted variation in fluorescence of Figure 6A, 
as a function of cycle number n. 

35 The following assumptions are preferably made: 



• the effectiveness of the reaction is relatively 
constant during the initial cycles; and 

after a certain number of cycles have been 
performed, the effectiveness of the reaction decreases. 

Figure 6B confirms the second assumption since it 
can be seen that effectiveness decreases as from cycle 
n=17. However, the measured effectivenesses in cycles 1 
to 16 are very noisy, which makes it difficult to verify 
the first assumption graphically. 

Nevertheless, it is preferable to assume that 
variation in effectiveness obeys a model of the type 
including : 

• a first stage that is constant between the first 
PCR cycle and the cycle (C EEP -1) preceding the changeover 
cycle written C EEP ;-and 

• a second stage in which it decreases for cycles of 
numbers greater than or equal to cycle (C EEP -1) . 

The cycle (C EEP -1) thus represents the last cycle 
(which may be a fraction) for which effectiveness 
continues to be constant. 

It is then proposed to model the effectiveness of 
the reaction as follows: 

E = |E 0 for 0 < n < (C EEP -1) 

n [ ( i + E o r p( - p(n - c - +1) - 1 for (C EEP -1) < n < (N-l) 

where Eq and P are real parameters which are estimated 
using the amplification curve of Figure 6A, or using the 
effectiveness curve of Figure 6B in a manner described 
below . 

In a variant, some other selection may be preferred, 
e.g. from the models Fl to F3 given below, particularly 
depending on the type of nucleic acid that is to be 
quantified. 



Fl : 


E n 


= exp(-P(n 


-Ceep+D ) -1 


F2 : 


En 


= exp(-p(n 


-C EEP +l) a ) -1 


F3: 


E n 


= a - exp ( 


-H(n-C nEP +l) a ) -1 



Preferably, several sets of parameters are estimated 
in step c) for several candidate changeover cycles C EEP , 
and the minimum candidate cycle is selected for which the 
associated parameters maximize the statistical 
correlations that can be undertaken in step c) , for each 
changeover cycle C EEP . 

As mentioned above, expression (1) may also be 
written in the form: 



n-l 

F n = Fo nd + E k ) for n e {l,2,...,N} (4) 

k=0 



Thus, by introducing the expression (3) for 
effectiveness into formula (4), a new model is obtained 
having four parameters (F 0 , E 0 , p, C EEP ) for the adjusted 



emitted fluorescence F n : 



F 0 (l + E 0 ) n for 1 < n < C EEP 

r ■ l-exp(-ft (n-C EEP ) ) ( 5 ] 

c EEP + 



F 0 (l + E 0 ) exp(&H for C EEP < n < N 



The initial size N 0 of the population of interest, 
the effectiveness E 0 of the reaction of n=0, the parameter 
P, and the changeover cycle C EEP are evaluated 
repetitively for several cycle values in the vicinity of 
the changeover region CHO in order to find a statistical 
correlation maximum that is achieved for a minimum cycle 
value that is equal to the changeover cycle C EEP . 

In this second implementation, it is preferred to 
model variation in the measured and adjusted quantities 
of fluorescence as a function of cycle number on the 
basis of the models or variation in effectiveness, and 
subsequently to carry out the correlations directly on 
the measured and adjusted quantities of fluorescence. 

It should be observed that by adjusting the measured 
emitted fluorescence for background noise, an artificial 
adjustment is also made on the initial fluorescence F 0 . 
Thus, estimating the parameters of the effectiveness 
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model on the basis of effectiveness measurements that are 
deduced from adjusted fluorescence measurements 
constitutes an additional source of error and it might be 
preferable to proceed in two stages as described below 
for the third implementation. 

Nevertheless, the second implementation as described 
is simpler and adapts well to PCR quantification using 
fluorescence measurements. it is based on the real 
measurements of fluorescence F f n which correspond to the 
fluorescence measurements adjusted for drift in 
background noise together with compensation e on said 
measurements. Once the background noise has been 

subtracted, we have a relationship of the following type: 



where s is a quantity that may or may not depend on cycle 
number n. It is preferably selected to be constant. 

Under such circumstances, the measured and 
"adjusted" effectiveness also written E' n on cycle n is 
defined by: 

E "=TTf - 1 = ~~f - 1 f °r all n e {1,2 N-l} (7) 

The model of above relationship (5) thus becomes: 



F' 0 (l + E 0 ) n for 1 < n < C E 

-exp(-p '(n-C EEP ) ) (8 ) 

EEP 



c .±1 

, , EEP 

(F 0 (l + E 0 ) 



'EEP 



exp(p,) - 1 for C FFP < n < N 



Under such circumstances, the effectiveness values 
E' n are approximated experimentally from the measurements 
so as to be able to set a minimum acceptable 
effectiveness threshold during the stage of decreasing 
effectiveness. A threshold cycle is thus determined 
beyond which the adjusted fluorescence measurements are 
not used for the purposes of the model (points NEG in 
Figure 4B) . Typically, the threshold cycle corresponds 
to the first cycle in the stage of decreasing 



effectiveness at which effectiveness drops below some 
minimum acceptable effectiveness threshold (e.g. 0.1E 0 ). 

More generally, the value of the effectiveness 
threshold preferably lies in the range 0 to 0.5, and PCR 
having an effectiveness value below said threshold is 
potentially biased by uncontrolled inhibition phenomena. 

In the example shown in Figure 7, the threshold 
value for effectiveness was set at 0.02 (i.e. 2% of E 0 ) . 
The threshold cycle C s corresponded to cycle n=36. 
Figure 7 shows the adjusted measurements of emitted 
fluorescence. It can be seen that there is satisfactory 
correlation with the model (continuous line) for those 
experimental measurements (marked with an "o") up to 
cycle n=36. Figure 8 also shows good correlation with 
experimental measu-rement s for predictive effectiveness as 
obtained from Figure 7 using the model based on measured 
and adjusted fluorescence. 

The main steps of this implementation can be 
summarized as follows, with reference to Figure 9. 

In a start step 70, the measured values for 
quantities of fluorescence have been obtained and 
adjusted relative to background noise as a function of 
cycle number n, as shown in Figure 6A. 

In step 71, an approximation for effectiveness of 
the reaction in cycle n is calculated using above formula 
(2) for each of the cycles n = 1, 2, (N-l). 

In step 72, the minimum cycle C s is determined for 
which the following two conditions are satisfied: 

the cycle C s lies in the stage of decreasing 
effectiveness ; and 

• the effectiveness of the threshold cycle is less 
than the threshold effectiveness value E s (e.g. E s = 
0.1E 0 ) : 
Ecs =S E s 

It is already possible to eliminate the points NEG 
for which effectiveness is less than E . 



In step 73, a model is formed for the curve of 
adjusted emitted fluorescence which effectiveness is 
decreasing over the cycle range C EEP = (C s - 5) to C s , 
using expression (8) in which it is assumed that 
compensation e is given by s = F' 0 : 

l-ex p(-ft(n-C EEP ) ) 
^EEP 

F n = F^l + E'J exp(PH1 - F' 0 

Thereafter, test 74 on the value E' 0 estimated for 
the value E f 0 and the decrementation in step 75 of the 
value for the changeover cycle C EEP seeks to find the 
looked-for value of C EEP using a step size P (which may be 
equal to 1) , and in repeating step 73 so long as the 
value of E' 0 is less than 1. 

Thereafter, when the estimated effectiveness value 
exceeds the value 1 (arrow n on exiting the test 74), the 
value of the index C EEP is incremented by a step of size h 
(which may be a fraction smaller than unity) in step 76 
and in step 77 fluorescence F n is modeled in the same 
manner as in step 73. So long as the estimated 

effectiveness E' 0 is greater than or equal to 1 in step 
78, steps 76 to 78 are repeated. When the estimated 
effectiveness takes a value of less than 1, the estimated 

/\ /v /v ^ 

parameters (F f 0 ,E f 0 ,p f 0 ,C EEP ) are conserved in an end step 
79. 

In this step, a value F 0 has finally been obtained 
that alone is representative of the initial population 
size N 0 in the sample of interest. It is then possible to 
use at least one standard sample having a known 
population size N 0 st so as to determine in step 80 the 
initial population size N 0 in the- sample of interest. 

For this purpose, a measured value of an initial 
population size F 0st in a standard sample of known initial 
population N 0st is obtained. Thereafter, the value of the 
initial population size N 0 in the sample of interest is 
obtained by deriving a proportionality relationship 
between the measurement for the standard sample and its 
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known initial population size, and applying that 
relationship to the measurement F' 0 to obtain the actual 
initial population size N 0 . 

In other words, in step 80 of Figure 9, it is 
5 possible to determine the value N 0 of the initial 
population size in the sample of interest by applying a 
simple proportionality relationship of the type: 
N 0 = F 0 (N Ost /F' 0st ) 

implying that the initial population size in the standard 

10 N ost and the ratio of the corrected fluorescences as 
compensated and estimated by adjusting the fluorescence 
model apply both to the sample of interest and to the 
standard sample. 

It will thus be understood that a single standard 

15 ought to be sufficient for determining the initial size 
of the population of interest in the sample of interest, 
which is an advantage provided by the invention. 

Nevertheless, in a variant, and where necessary, 
provision could also be made to obtain respective 

20 measured values for initial population sizes F' 0st in a 
plurality of standard samples having known initial 
population sizes N 0st . Thereafter, a dependency 

relationship is established between the initial 
population sizes N 0st of the standard samples and the 

25 respective measured values for their initial population 
sizes F' 0st . Thereafter, after finding the measured value 
for the initial population size of the sample of interest 
F' 0 , the actual initial population size N 0 of interest is 
determined by interpolation using the dependency 

30 relationship. It will be understood that this dependency 
relationship may also typically be a regression of the 
type shown in Figure 5, but having the initial 
fluorescence values F' 0st and F' 0 of the standards and of 
the sample of interest plotted up the ordinate (or the 

35 values of their respective logarithms) instead of 
plotting values for the changeover index C EEP . 
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Once use is made of one or more standards, provision 
can be made for one or more standard samples having 
respective known initial population sizes N 0st to which 
the succession of amplification reactions is applied 
5 under substantially the same conditions as for the sample 
of interest. Thereafter, the measured values F 0st for 
their initial population sizes are determined by making 
comparisons of the parameterized variations with the 
experimental measurements, as for the sample of interest. 

10 In other words, the same calculations are naturally 

applied concerning the measured and adjusted quantities 
of fluorescence both on the standard (s) and on the sample 
of interest. The quantity of fluorescence F 0st before the 
first cycle is estimated for the standard (s) using the 

15 same method as is -used for determining F 0 for the sample 
of interest, as described above. 

A third implementation, corresponding to a variant 
of the above-described second implementation consists 
overall in adjusting the model for the effectiveness E n 

20 relative to the experimental measurements, and in 
subsequently injecting said adjusted effectiveness model 
into the model for the current population size N n , or into 
the model for the measurement F n . This third 

implementation can be summarized as follows. 

25 The parameterized variation constructed in step b) 

is representative of yield, and in step c) , experimental 
variation of the yield is determined on the basis of 
experimental measurements in order to compare the 
parameterized variations with the experimental variation. 

30 Thereafter, in order to obtain a parameter representative 
of the initial population size N 0 the following steps are 
performed in step d) : 

dl) determining a second parameterized variation 
representative of the current population size N n in the 

35 sample of interest, making use at least of the parameter 
representing the changeover index C EEP , and a parameter 
representative of the initial population size N 0 ; r 
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d2) applying to said second variation, a 
parameterized value for the changeover index C EEP as 
determined in step c) ; and 

d3) adjusting at least the parameter representative 
5 of the initial population size N 0 by direct comparison of 
the second variation with the experimental measurements. 
Advantageously the following are performed: 

in step d2), applying a coarse value for the 
changeover index C EEP in the same manner as described for 
10 detecting it with reference to above Figure 4B; and 

• in step d3), subsequently refining the index value 
together with adjusting the parameter representative of 
the initial population size N 0 . 

Finally, it should be understood that the presently 
15 preferred second - implementation of Figures 7 and 8 
differs from this third implementation by the fact that 
no attempt is made to perform correlations on 
effectiveness, but use is made merely of the mathematical 
model for effectiveness variation in order to model and 
20 refine the estimate of the correct and compensated 
fluorescence . 

Naturally, the present invention is not limited to 
the embodiments described above by way of example, and it 
extends to other variants. 

25 Thus, it will be understood that the present 

invention can also apply to relative quantification, in 
particular by PCR. In this application, as well as 
amplifying the population of interest, a reference 
population is also amplified either simultaneously in the 

30 same medium, or separately. Measurements are taken as 
follows : 

experimental measurements representative of the 
size of the population of interest; and 

experimental measurements representative of the 
35 reference population size. 

The method can then continue by applying steps a) , 
b) , and c) to the reference population while step d) 
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consists merely in determining a ratio between the 
respective initial sizes of the population of interest 
and of the reference population. 

Relative quantification can be used for analyzing 
the expression of a gene of interest during the 
development of an organism. In order to correct in 
particular for variations in quantity and in quality 
between samples taken from the organism at different 
times, in addition to analyzing the target gene of 
interest, a reference gene is also analyzed that is known 
for having a level of expression that remains stable 
during development . 

A final step then consists in comparing the ratios 

^Otarget 



N 0ref 



between the various samples that have been taken. 

In order to achieve the desired results, two 
strategies are possible. 

The prior art strategy is based on detecting the 
threshold cycle Ct and it normally takes place as 
follows. For each sample taken at different instants tO, 
tl, t2, tn, the ratio 

^Otarget 



N Oref 



is determined, making use of at least one standard (i.e. 
a standard for which N otarget and N 0ref are known) , which 
amounts to performing two successive absolute 
quantifications followed by calculating a ratio. 

Another strategy that is particularly advantageous 
in the context of the invention consists in determining 
for each sample taken at different instants tO, tl, t2, 
. . . , tn the ratio : 

( ^Otarget 
^Oref y sample 



^Otarget 
fc ^Oref / S ample_tO 



directly by using the following formula: 



Otarget 



Otarget 'sample 




N, 



Oref 



Otarget ^sample_tO 



sample 



Otarget 



Oref 'sample 



V Oref /sample_tO 



Oref 'sample_tO 



In this second implementation, which in the end 
makes use only of the parameter F 0 , in combination with 
the technique of the invention, no standard sample is 
needed, which is particularly advantageous. 

Reference is now made to Figure 10 which shows an 
installation for implementing the method of the 
invention. It comprises a support SUPP in this case 
comprising a well ^ containing the sample of interest ECH 
and a well containing a standard sample referenced St, 
for example. The support SUPP is enclosed in an 

enclosure ENC, e.g. fitted with heater means (not shown) 
for applying a PCR reaction to the standard and to the 
sample of interest. 

In the example described, provision is preferably 
made to take measurements of the quantities of 
fluorescence emitted on each cycle, both by the standard 
St and by the sample of interest ECH. To this end, a 
selected reagent is inserted into the wells and the 
samples are illuminated by a lamp (e.g. a halogen- 
tungsten lamp) in order to measure the respective 
quantities of fluorescence coming from the sample of 
interest and from the standard sample on each PCR cycle 
that is applied. In addition, an apparatus for detecting 
fluorescence comprises, for example, an objective lens 11 
for collecting the light coming from the fluorescence, 
and photon counting means 10, e.g. a charge-coupled 
device (CCD) camera, and/or photomult iplier s , in order to 
measure the fluorescence emitted on each PCR cycle from 
the sample of interest and from the standard. Thus, the 
fluorescence emitted by each well is advantageously 
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focused by the lens 11 and then is preferably detected by 
a CCD camera 10 connected to an acquisition card 21, e.g. 
of the Personal Computer Memory Card International 
Association (PCMCIA) type provided in a central unit 20 
5 of a computer. 

The computer is then connected to the above- 
mentioned measuring means 10 to receive therefrom signals 
that are representative of the measured quantities of 
fluorescence detected on each PCR cycle, and to process 
10 these signals in order to determine an initial size for 
the population of interest prior to the first cycle, by 
implementing the method of the invention. 

Typically, the processor unit comprises the 
following : 

15 • an acquisition card 21 connected to the 

measurement means 10; 

working memory 25 (e.g. of the random access 
memory (RAM) type) for temporary storage and processing 
of the above-mentioned signals; 

20 • permanent memory 24 for storing the computer 

program product in the meaning of the invention and for 
storing the data that has been processed and that is 
ready for use, e.g. in subsequent diagnosis; 

• where appropriate, a reader 22 of a memory medium 
25 such as a CD-ROM, for example, which may initially carry 

the computer program product; 

optionally a communications interface 26 for 
communicating with a local or remote site (connection 
28), e.g. for transmitting the processed data so as to 
30 enable a diagnosis to be made remotely concerning a 
patient ; 

• a graphics interface 27 typically connected to a 
display screen 30; and 

a processor 23 for managing the. interactions 
35 between these various items of equipment. 



The computer may also have input members such as a 
keyboard 41 and/or a mouse 42 connected to the central 
unit 20. 

Nevertheless, it should be understood that in the 
meaning of the invention the installation comprises 
overall : 

• a sample support SUPP, at least for the sample of 
interest ; 

• a first apparatus ENC for applying said succession 
of amplification reactions at least to the population of 
interest in the sample of interest; 

a second apparatus 10 for taking measurements 
representative of the current size of the population of 
interest; and 

computer . means 20 suitable for receiving 
measurement signals from the second apparatus 10 and for 
implementing all or some of the steps of the method of 
the invention. 

For this purpose, a computer program product can be 
used for controlling the computer means. The program may 
be stored in a memory of the processor unit 20 or on a 
removable memory medium (CD-ROM etc.) and suitable for 
co-operating with the reader of the processor unit. The 
computer program in the meaning of the invention then 
contains instructions for implementing all or some of the 
steps of the method of the invention. For example, the 
algorithm of the program may be represented by a flow 
chart equivalent to the diagram of Figure 9. 



